y = [2 2; 3 2; 5 6; 2 8; 9 2; 11 12];
barHdl = bar(y);
CList = [87,87,213; 138,213,95]./255;
for i = 1:length(barHdl)
    barHdl(i).FaceColor = CList(i,:);
end

legend()
set(gca, 'XLim',[.5,6.5])
SAxes()


xregion(.5,3.5, 'FaceColor',[233,241,254]./255)
xregion(3.5,6.5, 'FaceColor',[251,244,218]./255)

bar2_5D(barHdl, 'k')
